Determining orbits for the Milky Way's dwarfs 
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Abstract. We calculate orbits for the Milky Way dwarf galaxies with proper motions, and compare 
these to subhalo orbits in a high resolution cosmological simulation. We use this same simulation to 
assess how well are able to recover orbits in the face of measurement errors, a time varying triaxial 
gravitational potential, and satellite-satellite interactions. We find that, for present measurement 
uncertainties, we are able to recover the apocentre r a and pericentre r p to ~ 40%. However, even 
with better data the non-sphericity of the potential and satellite interactions during group infall make 
the orbital recovery more challenging. Dynamical friction, satellite mass loss and the mass evolution 
of the main halo play a more minor role. 

We apply our technique to nine Milky Way dwarfs with observed proper motions. We show that 
their mean apocentre is consistent with the most massive subhalos that form before z = 10, lending 
support to the idea that the Milky Way dwarfs formed before reionisation. 
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Testing the method. In [1] we derive orbits for the Milky Way's (MW) dwarfs with 
proper motion data, by integrating their orbits backwards in a fixed potential. We test 
this procedure using the high resolution Via Lactea I cosmological simulation of a MW 
analogue [VL1; 2]. For this we extract from the simulation two sets of subhalos: the 50 
most massive today (Zq ), and the 50 most massive before redshift z = 10 (z^q) 1 . In both 
cases, we include only subhalos with mass M > 1O 7 M and distance to the centre of the 
main halo r < 150 kpc at redshift z = 0. 

To test the limitations of this approach, we integrate the orbits in 3 different models: 
The fiducial model: This uses a static, spherical, NFW potential [3] fit to the VL1 main 
halo at z = [2]; The dynamical friction model : This uses the same potential with 
dynamical friction forces [4], a prescription for the mass loss of satellites as well as 
the mass growth of the main halo [5]; The triaxial model : This uses a triaxial NFW 
potential as in [6, 7]. 

We find that the orbits without measurement errors are systematically effected by an 
incorrect halo shape as well as satellite- satellite interactions. Dynamical friction, mass 
loss of the satellite and mass evolution of the main halo play a minor role. Current 
measurement errors are more significant than model systematics. They bias the mean 
apo-/pericentre to higher values. We are able to recover apo-/pericentre up to 40%. 



See the data at http://www.ucolick.org/diemand/vl/ 
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FIGURE 1. Recovered peri-/apocentres for nine Milky Way dwarfs with observed proper motions. The 
error bars show the values for the oblate potential of [8]. Overlaid are the mean and standard deviation of 
the values from Zq° (light grey band, dashed line) and z\q (dark grey band, solid line). The black diamonds 
denote the values derived from the Sagittarius stream [8]. 



Milky Way Dwarfs. We derive orbits for the MW dwarfs with proper motion data, 
by integrating their motion backwards in the oblate potential from [8]. We take current 
measurement errors into account by building an ensemble of 1000 orbits for each dwarf 
drawn from its error distributions. The results are shown in Figure 1. Our recovered 
mean pericentre (r p ) is higher than the mean pericentre in the simulation. This could 
be either due to satellite depletion by the galactic disk as recently discussed by [9]; but 
can also be explained by the bias we see in the VL1 recovered data which is due the 
large proper motion errors. This can be reduced by better quality data. The apocentre 
distances r a instead are lower than the mean of the zjj sample, but consistent with the 
z\q sample. This lends further support to the idea that the MW's dwarfs formed early 
before reionisation [e.g. 10-12]. 



REFERENCES 

1 . H. Lux, J. I. Read, and G. Lake, in prep. (2010). 

2. J. Diemand,M. Kuhlen, andP. Madau, ApJ 657, 262-270 (2007), arXiv : astro-ph/0 61137 0. 

3. J. F. Navarro, C. S. Frenk, and S. D. M. White, ApJ 462, 563-+ (1996), arXiv: astro-ph/ 
9508025. 

4. S. Chandrasekhar, ApJ 97, 255-+ (1943). 

5 . H. Zhao, MNRAS 351, 89 1-902 (2004), arXiv:astro-ph/0403564. 

6. M. Kuhlen, J. Diemand, and P. Madau, ApJ 671, 1 1 35-1 146 (2007), 0705.2037. 

7. J. Guedes, P. Madau, M. Kuhlen, J. Diemand, and M. Zemp, ApJ 702, 890-900 (2009), 907 .08 92. 

8. D. R. Law, K. V. Johnston, and S. R. Majewski, ApJ 619, 807-823 (2005), arXiv: astro-ph/ 
0407566. 

9. E. D'Onghia, V. Springel, L. Hernquist, and D. Keres, ArXiv e-prints (2009), 0907.3482. 

10. A. V. Kravtsov, O. Y. Gnedin, and A. A. Klypin, ApJ 609, 482^97 (2004), arXiv: astro-ph/ 
0401088. 

11. B. Moore, J. Diemand, P. Madau, M. Zemp, and J. Stadel, MNRAS 368, 563-570 (2006), arXiv : 
astro-ph/0510370. 

12. A. V. Maccio', X. Kang, F. Fontanot, R. S. Somerville, S. E. Koposov, and P. Monaco, ArXiv e-prints 
(2009), 0903 . 4681. 



